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ABSTRACT 

We present results from an extensive survey of 64 cavities in the X-ray hales of clusters, groups 
and normal elliptical galaxies. We show that the evolution of the size of the cavities as they rise 
in the X-ray atmosphere is inconsistent with the standard model of adiabatic expansion of purely 
hydrodynamic models. We also note that the majority of the observed bubbles should have already 
been shredded apart by Rayleigh- Taylor and Richtmyer-Meshkov instabilities if they were of purely 
hydrodynamic nature. Instead we find that the data agrees much better with a model where the 
cavities are magnetically do minated a nd inflate d by a current-dominat ed magneto-hydrodynamic jet 
model, recently developed bv lLi et al] (j2006. ) and lNakamura et al] (|2006D . We conduct complex Monte- 
Carlo simulations of the cavity detection process including incompleteness effects to reproduce the 
cavity sample's characteristics. We find that the current-dominated model agrees within Icr, whereas 
the other models can be excluded at > 5cr confidence. To bring hydrodynamic models into better 
agreement, cavities would have to be continuously inflated. However, these assessments are dependent 
on our correct understanding of the detectability of cavities in X-ray atmospheres, and will await 
confirmation when automated cavity detection tools become available in the future. Our results have 
considerable impact on the energy budget associated with active galactic nucleus feedback. 

Subject headings: cooling flows — X-rays: ISM — X-rays: galaxies: clusters — galaxies: clusters: 
general — instabilities — MHD 



1. INTRODUCTION 

One of the current most active areas of research in 
X-ray astronomy is aimed at solving the problem of 
feedback in elliptical galaxies, groups and clusters of 
galaxies. All of these systems suffer from the same un- 
solved problem that the gas in their cores has cooling 
times significantly shorter than their age, suggesting that 
the gas should already hav e cooled and formed stars 
(|McNamara fc NulsenI (20071 ). This problem, historically 
referred to as the "cooling flow problem" , points toward 
the existence of a heat source that is preventing the gas 
from cooling effectively in the center. Due to a lack of 
better understanding, this effect is generally referred to 
with the generic term "feedback" . 

Many possible feedback mechanisms have been pro- 
posed in recent years. The most basic idea is the conver- 
sion of gravitational energy into thermal energy, as gas 
falls into the center towards a deeper gravitational poten- 
tial. However, this mechanism has been shown to only be 
significant in elliptical galaxies with a rather steep gravi- 
tational potential, such as NGC 6482 (jKhosroshahi et ahl 
[2004). and does not solve the cooling flow problem. 

Another heat source that has often been suggested 
is feedback from supernovae. If the gas cools suffi- 
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ciently, stars will form, which will then produce su- 
pernovae, which in tur n will heat the ambient medium 
(|Binnev fc Taboil ITOOl ). While this effect is certainly 
important, the amount of energy released can only be 
sufficient to offset cooling in low-luminosity elliptical 
galaxies, but does not even come close to the radiative 
losses in groups or cluster, or even high-luminosity el- 
lipticals. Even in normal elliptical galaxies, supernova 
feedback models have been shown to need fine-tuning 
to produce gas halos that r esemble X-ray observations 
(|Mathews fc Brighentill2003D . 

Instead, the central active galactic nucleus (AGN) has 
emerged as the most popular heat source in the last few 
years, though this still remains an active ma tter of de- 
bate in the community. iMcNamara fc NulsenI |2GG7) give 
a recent review of AGN heating i n clusters. Chandra and 
XMM-Newton observations (e.g. iMcNamara et al.|[2000f) 
have revealed the existence of giant cavities in clusters, 
which are often associated with radio emission from on- 
going jet activity of the central AGN, actively inflating 
these "bubbles" . This dir ect link has also be en observed 
on the small er group (e.g.lMorita et al.ll2006l) and galaxy 
scales (e.g. jFinoguenov fc Jones! l2001t iDiehl fc Statled 
Hoc?, 2008a3I3). AGN outbursts have been shown to have 
potentially more than sufficien t energy to offset c ooling 
in a large fraction of clusters (jBirzan et all 120041 ) ■ and 
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even in nor mal elliptical galaxies and groups ()Best et al.l 
[2005l[2?)06h . 

However, while the AGN has potentially sufficient en- 
ergy to offset cooling, the detailed dynamics of how the 
released energy can actually be transferred into ther- 
mal energy of the gas is poorly understood. A lot of 
theoretical work has been undertaken in the purely hy- 
drodynamical regime (e.g. lBriiggenll2003D to understand 
jet dynamics in inflati ng these cavities, but with h mited 
success so far (see e.g. IVernaleo fc Revnoldsl [20061 ). The 
most realistic simu lations to date have been conducted by 
I Heinz et"al] (|2006( l. who employ the dentist's drill model 
combined with a dynamical cluster initial setup to avoid 
some of the issues other more idealized simulations have. 
But it remains to be seen how successful this model is 
on the less dynamic group and galaxy scale, where AGN 
heat ing has recently been found to be important as well 
/e.g.lDiehl fc St atler''2007'. '2008a'; 'B est et al .1120051 [20061 : 
lAllen et al.l l2006: Shurkin ct al. 2007^ 

A problem of purely hydrodynamic jet simulations is 
that they ubiquitously produce a shock at the termina- 
tion point of the jet, which results in hot, shocked gas 
surrounding the jet lobe. This is contrary to the vast 
majority of X-ray observations, though shocks around 
radio lob es ha v e rece ntly been reported in NG C 3801 
dCros t on et al.l l2007f ) and Centaurus A (Kra ft et all 
l2003l l2007f). Another major issue is that the cre- 



ated cavities often do not morphologically resemble 
X-ray observations, and get shredded apart quickly. 
Rayleigh- Taylor and Richtmyer-Meshkov instabilities are 
expected to grow on the top of the bubble and Kelvin- 
Helmholtz instabilities develop on the shear flow at 
the bubble sides, as it rises. Nevertheless, bubbles 
with rise times on the order of several times the pre- 
dicted instability timescales have been observed, call- 
ing this scenario into question. Something has to sup- 
press the growth of th ese instabilities. In gene ral, tan- 
gen tial magnetic fields llJones fc De Young|[2005l) . viscos- 
ity (IRevnolds et al.l 12005), or continuous inflation of the 
bubble (IPizzo lato et al.ll20d^ have been proposed to sup- 
press the growth of instab il ities. 

Alternativ ely, iLi et al.l (|2006D and iNakamura et al.l 
(|2006l l2007f l have proposed a magnetically dominated 
jet model instead. Their 3D magneto-hydrodynamic 
(MHD) simulations advocate a current-carrying jet that 
evolves purely by injecting magnetic flux and energy 
into a small volume surrounding the central black hole. 
The injected magnetic fields are not force-free and inter- 
act with the ambient medium to form a self-coUimated 
jet. As the background pressure drops, the jet expands 
into a wide jet lobe, excavating a bubble in the ambi- 
ent medium. Mock X-ray observations of these systems 
show a s trong resemblance to observations of cavities in 
clusters (|Diehl et al.ll2007f ). In particular, these cavities 
generally expand subsonically and appear cooler than the 
surrounding medium, similar to observations. 

This paper addresses the question of whether one can 
use the bubbles' sizes as they rise in the intracluster 
medium as a function of their distance to the cluster 
center as an indicator for which model may indeed be 
correct. 

The paper is organized as follows: In section [2] we de- 
scribe the methodology on how we identify the cavities, 
and extract their main properties. fJSjthcn presents the 



predictions of the various models that we are consider- 
ing, and focus on their predictions for bubble sizes as a 
function of radius. In section 3) we show how this relates 
to observations. We then discuss various incompleteness 
issues associated with the data in sectionjU before we try 
to reproduce the properties of our data set with Monte- 
Carlo simulations in SjHland conclude (SJ7]). 

2. X-RAY DATA 

2.1. Cavity Sample 

Our cavity sample is mostly congruent with that of 
iRaffertv et al.l 12006). The only differences are that we 
exclude Abell 1068, which does not contain a ny cavities, 
and ad d the multiple cavities in Hydra A from lWise et al.l 
(|2007t l. Overall, the sample comprises a set of 64 individ- 
ual cavities in 32 different hosts. Most of the hosts are 
clusters (30/32), with one group (HCG 62) and one in- 
dividual galaxy (M 84) . The cluster sample incorporates 
very poor clusters (e.g. Abell 262) to very rich clusters 
(e.g. Abell 1795). For clarity and simplicity, we will refer 
to the hosts of our cavity sample in general as "clusters" , 
despite the two exceptions noted above. All of conclu- 
sions are not changed even if we exclude the one galaxy 
and the one group. 

2.2. Data Analysis 

The X-ray data, used to derive the cavity properties, 
are taken from the Chandra Dat a Archive and repro- 
cessed in CIAO 3.3, as described in lRaffertv et al.l ()2006f) . 
In summary, we start from level 1 event files and apply 
a uniform calibration across the sample using CALDB 
3.2.0. All systems have been observed with the ACIS de- 
tector in non-grating mode. During reprocessing, correc- 
tions to account for the time-dependent charge transfer 
inefficiency (CTI) and gain problems are applied. Ad- 
ditionally, the data are filtered for bad grades and hot 
pixels, and cosmic-ray events are identified and flagged. 
Lastly, exposure maps are made using weights appropri- 
ate for a single-temperature plasma. 

Blank-sky background files^ are used to estimate the 
background, as the cluster emission often fills the entire 
field of view of the CCD. Periods of high background 
are excluded from the source observation by filtering the 
data for periods during which the count rate exceeds 20% 
the quiescent background count rate, after exclusion of 
all sources. The blank-sky background rate is normal- 
ized by the count rate in the 10-12 keV energy range and 
reprojected to match the aspect solution of the source 
data. After reprojection, the background spectra are ex- 
tracted from the same regions of the CCD as the source 
spectra. 

2.3. Extraction of Cluster Properties 

For this work, we additionally extract X-ray surface 
brightness profiles using elliptical annuli centered on the 
X-ray core. The geometry of the annuli is fixed to that 
of the best-fit annulus at an intermediate radius, gener- 
ally beyond the location of the cavities or other structure 
where the cluster emission is fairly smooth. The profiles 
are extracted between 0.1-12 keV and corrected for ex- 
posure with the exposure maps. 

^ See |http: / / cxc.harvard.edu/contrib/maxim/bg/ 1 
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The surface brightness profiles are then logarithmically 
binned to avoid that the model fits are driven mainly by 
noisy data at large radii. We fit a /9-model profile with an 
additional background constant to each surface bright- 
ness profile by minimizing the x^-value with the publicly 
available IDL^ fitting package MPFIT ^. For more infor- 
mation on the /3-model, please refer to ISarazinI (|1988f) or 
section [3Tl Host properties are summarized in Table [TJ 

2.4. Extraction of Bubble Properties 

We identify and measure 58 bubbles in the X-ray data 
(in addition to th e the six cavities in Hydra A identified 
and measured bv lWise et al.ll2007h . Each bubble is con- 
firmed as likely to be associated with the central AGN 
by compa ring its location t o that of the AGN's radio 
lobes (see lBirzan et al.ll2008[ ). We note that most of the 
bu bbles have beeii previo usly identified in the literature; 
see iRafFertv et al.l ()2006[ ) for references. In most cases, 
there is a direct anti-correlation between the X-ray and 
radio emission, such that the radio lobes fill cavities in 
the X-ray emission. 

We measu re the bubble p roper ti es following the proce- 
dure used in Birzan et~all (|2CI04D : iRaffertv et al.l (|2006D . 
Briefly, we estimate the size and position of each cavity 
by overlaying ellipses onto the exposure-corrected, un- 
smoothed images. The ellipses are positioned such that 
they encompass the cavity, but not any bright emission 
associated with the cavity rim, if present. 

To obtain the density and temperature at the bub- 
ble's projected location, we perform a deprojection of 
the cluster emission using the PROJCT model in XSpec 
11.3.1. Spectra of at least 3000 counts in the 0.5-7 
keV band are extracted in concentric annuli, along with 
the associated weighted responses, using standard CIAO 
tools. The spectra are then fit with models of a single- 
temperature plasma (MEKAL) with foreground absorp- 
tion (WABS) between the energies of 0.5-7 keV. The 
foreground co lumn density is held fixed to the Galac- 
tic value of ,Dickev fc Lockma n (1990), and the redshift 
of the MEKAL comp onent is fixed to the values given in 
IRaffertv et al.l ()2006f ). The tei nperature and a bundance 
(relative to the solar values of I Anders fc Grev esse 1989) 
of the MEKAL component are allowed to vary freely in 
each shell, and the density is calculated from the normal- 
ization of the MEKAL component assuming = 1.2nH. 
Finally, the pressure is calculated as p = nkT, where we 
assume a fully ionized ideal gas. The ambient pressure 
and temperature at the bubble's location are derived by 
interpolating the pressure and temperature profiles at the 
projected radius of the bubble's center. Bubble proper- 
ties are summarized in Table [5] 

3. THEORY: PREDICTIONS FOR BUBBLE SIZES 
AS A FUNCTION OF RADIUS 

3.1. Assumptions and Underlying Thermodynamic 
Structure 

To distinguish between different models of X-ray cavi- 
ties, we need to compare definite predictions of the vari- 
ous models. One such well-defined prediction is the evo- 
lution of the bubble sizes as they rise in the intracluster 



medium. Since X-ray observations show no evidence for 
shocks around cavities, we assume throughout this pa- 
per that their inflation and expansion is subsonic and in 
approximate pressure equilibrium with the surrounding 
gas. In the following sections we will summarize the pre- 
dictions of various models and show their predictions for 
the bubble expansion as a function of radius. 

To make these predictions, we have to assume a back- 
ground pressure profile, which we choose to be the widely 
known and used /3 model fe.g. lSarazinlll988l ). This model 
adequately describes a wide range of cluster profiles over 
a rather wide range of radii, and its shape is uniquely de- 
fined by only two parameters: the core radius Tc and the 
exponential parameter /3, which determines the sharp- 
ness of the turnover beyond the core radius, as well as 
the asymptotic slope of the profile. The surface bright- 
ness profile can be described as 

S{R) = S^[l + {Rlr,)Y^^-'\ (1) 

where R denotes the projected radius."' 

One of the most advantageous properties of the (3 
model is that one can analytically derive its deprojec- 
tion into three dimensions for the case of perfect spher- 
ical symmetry, which has made this model the de-facto 
standard model used for clusters in X-ray astronomy. As- 
suming isothermality at a temperature T throughout the 
cluster then allows us to determine the density and pres- 
sure profile as a function of the spherical radius r (e.g. 
lEttorill200l : 

_ 3/3 

p(r) = po [1 + {rlr,f] ^ , and (2) 

_ 3/3 

Pgas(r)-po[l + (r/rc)2] (3) 

Here, the central pressure po is connected to the central 
density pg and isothermal temperature T and average 
molecular weight fh via the equation of state of an ideal 
gas with an adiabatic index F = 5/3: po = Po'fn~^ kT. 
Though this represents a simplified picture of a real clus- 
ter, it allows us to make analytical predictions which 
should at least qualitatively hold up even in more com- 
plex systems. Real clusters tend to exhibit cool cores, an 
outwardly rising temperature gradient, and sometimes 
also an additional surface brightness excess near the cen- 
ter. We did repeat the analysis that follows with the 
more complex double-/? and temperature models, which 
leaves our conclusions unchanged and results in no ad- 
ditional insight into the problem. Thus, we decided to 
discuss only the isothermal single-/? model for simplicity. 
About 10 clusters in our sample show signs of a second 
(3 model in the very center. However, we made sure that 
the single-/? model fit adequately represented the pres- 
sure profile around the location of the cavity, which is all 
our discussion here requires. 

In the following, we will present predictions of various 
models with regard to the size evolution of cavities as 
they rise in this ambient pressure profile. We assume 
that they are always in approximate pressure equilib- 
rium with the surroundings. This assumption is justified 
by the general absence of strong shocks surrounding the 
bubbles, though there is now evidence for weak shocks 



^ http: //www. ittvis . com/idl/ 

^ http: //cow. physics .wise . edu/~craigni/idl/f itting.html 



^ T his corresponds to the profile assumed in INakamura et aO 
II2007I ') with K = 3/3/2 as used in their MHD simulations. 
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TABLE 1 
Cluster X-ray Properties. 
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with Mach nu mbers ranging from 1.2-1.7 in a number 
of systems (see iMcNamara fc NulsenI 120071 for a recent 
review) . These weak shocks are usually not directly asso- 
ciated with the cavity rims however, but may rather indi- 
cate an earlier faster expansion phase of the bubble. Only 
two objects are known with evidence for shocks in direct 
contact wi th the radio lob es: Centaurus A (K r aft et al.l 
[200112001 and NGC 3801 (|Croston et all^OOl . In fact, 
the majority of rims around bubbles seem to be cooler 
than the ambient medium, suggesting a subsonic infla- 
tion and rise of the bubbles, which provides a solid basis 
for our assumption. 

3.2. Adiabatic Expansion of Purely Hydrodynamic 
Bubbles (AD43/AD53) 

The standard assumption in most hydrodynamic sim- 
ulations of cluster-scale jets is that the jet is kinetically 
dominated and inflates a bubble at the location of the 
jet lobes. We assume that this bubble is filled with 
ideal gas with an adiabatic index F. For simplicity, we 
will assume that the bubble is inflated at the center for 
this model with an internal pressure pb^ and Volume 
Vb.o — 4/37ri?^Q. How the total outburst energy Er is 



related to the mechanical work associated with inflating 
this bubble, is a matter of debate and can be parame- 
terized by an efficiency parameter r/ via Er = Jyp&.o^^.o- 
Realistic 3-dimen sional hydrodynamic simulations (e.g. 
iHeinz et al.]l2006D are needed to determine the value of 77 
and to find out whether it is a function of environmental 
parameters, cluster dynamics, black hole mass, or similar 
quantities. The initial bubble radius Rtfi at the cluster 
center should be a direct function of the outburst energy 
and this efficiency parameter, such that 

«'-(^)"\ « 

Our analysis is insensitive to the details of this mech- 
anism. For now, let us simply assume that we have a 
purely hydrodynamic, inflated bubble which adiabati- 
cally expands while rising in the intracluster medium. 
This bubble should expand according to the well-known 
adiabatic law ptVj^ = pb^Vj^Q. Combining this relation 
with the assumed thermodynamic structure of the clus- 
ter gas predicts the evolution of the bubble radius Rb.r 
as a function of radius, assuming the bubble stays in- 
tact and is constantly in pressure equilibrium with the 



The Nature of X-ray Cavities 



TABLE 2 
Cavity X-ray Properties. 
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Note. — Column H-2 == Cavity number and host name; Column 3 - Projected distance of the bubble center 
to the centroid of the cluster X-ray emission; Column 4 - Size of the bubble; Column 5 - pV work associated with 
inflating a spherical bubble of size V = 4/37rij3 at r; Column 6 - Current for CDJ model (from equation llll l: 
Column 7 — Energy injection rate for the CIH model (from equation[6]l; Column 8 — Bubble age assuming a constant 
rise speed at the speed of sound; Column 9,10,11 - Kelvin-Helmholtz, Ray leigh- Taylor and Richtmyer-Meshkov 
instability timescales for the largest possible mode (Rt)- 
* Hydra A data taken fromlWise et all (120071) 
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ambient medium: 
Rb.r 



Rb^o [1 + (r/rc) 



2] ^ 



(5) 



Thus, at large radii, the cavities should grow asymptoti- 
cally as Rbx {r/rc)~^^^ ■ For typical values of P — 0.5 
for a cluster and a relativistic gas with F = 4/3, we get 
Rbx oc [r/r^f/^] for F = 5/3 we get i?b,r oc [r/r^f/^^. 
From here on, we will refer to this model as the adiabatic 
expansion model, or short AD43 and AD53 for F = 4/3 
and 5/3, respectively. Both correlations are actually sur- 
prisingly shallow, especially as one should rather con- 
sider these bubble size predictions as upper limits, since 
they do not take the effects of hydrodynamic instabilities 
into account. In fact, we will show later in section [5.31 
that these instabilities should have shredded most hy- 
drodynamic bubbles apart. Hydrodynamic simulations 
including viscosity have recently been found to stabilize 
the bubbles and suppr ess the development of instabilities 
(|Revnolds et al.ll2005D . 

Another effect not included in these predictions is the 
radiative loss of pressure support. If the highly relativis- 
tic particles inside the cavities provide pressure support 
and radiate a significant fraction of their energy away, 
the predicted size evolution of the cavities would be even 
shallower. Typical radiative lifetimes for the relativis- 
tic electrons are around 10^ — 10^ yr, o n par with bub- 
ble ages. However , it ha s been shown (|Dunn fc FabianI 
120041 : iBirzan et al] l2008l ) that the relativistic electrons 
are insufficient to supply the pressure needed to keep the 
bubbles in pressure equilibrium. The missing pressure 
component must be either supplied by "extra" particles, 
such as protons or entrained thermal gas, or magnetic 
fields. 

3.3. Continuously Inflated Hydrodynamic Bubbles 
(CIH) 

The just described AD43 and AD53 models both as- 
sume the existence of an a priori preformed bubble at 
the center which then buoyantly rises. In this extreme 
case, the outburst energy is deposited instantaneously 
by the AGN. In this section, we now consider the other 
extreme with a continuous energy injection over time. 
For simplicity, we assume this process to be constant 
(i.e. E{t) = Ei^t). Assuming the AGN injects purely 
hydrodynamic material, we can rewrite this as a change 
in initial volume Vq with efficiency 77 over time, rather 
than injected energy: E^t = r]poVo{t) = 4:/3r]poTrRl{t). 
Replacing Rb.o in equation [5] with its new time depen- 
dent form, and assuming a continuous rise speed (e.g. at 
the sound speed, Cg = r/t), yields the prediction for the 
radial evolution of the bubble size for the continuously 
inflated hydrodynamic (CIH) bubble model: 



R 



6. CIH 



y AtttjpqCs 



1/3 



[1 + (r/r,)2 



(6) 



Asymptotically at large radii, this model predicts that 
the bubble size scales as Rb cx r^^/^^ for a F = 4/3 and 
a (3 parameter of 0.5. 

Note that this model implicitly relies on several ba- 
sic assumptions. First, we assume in equation ([6]) that 



all of the new material that is injected by the AGN is 
directly incorporated into the cavity at radius r, and in- 
stantly adiabatically expanded to match the outside gas 
pressure. In a realistic model, this would require the 
jet "nozzle" to always be directly attached to the ris- 
ing bubble, while keeping the bubble intact. Realistic 
3D hydro simulations are needed to address the ques- 
tion whether this is actually feasible, without produc- 
ing a hot spot for the majority of cavities. Current 
state -of-the-art simulations tend to suggest otherwise 
(e.g. IVernaleo fc Revnoldsl [200I [BTliggcn 20Cg). An- 
other problem of this model is the length of the injec- 
tion time, which in our toy model would then be equal 
to the bubble ages trisc given in Tabled These ages are 
are about an order of magnitude longer than what is ex- 
pected for typical AGN lif e times and duty cycles (e.g. 
lAdelberger fc Steide]|[200l . 

Thus, this CIH model should rather be considered the 
other "extreme" and realistic paths in the r — Rb diagram 
should fall in between the CIH and AD models. In fact, 
one could realize arbitrary injection functions Ea{t), and 
accordingly replace the injection term in equation [6l 

3.4. Magnetically Dominated Bubbles with Flux-Frozen 
Magnetic Loops (FML) 

In the following, we show the prediction for the size 
evolution of a magnetically dominated bubble for the 
case of a flux-frozen magnetic dipole field, refering to 
this model as the FML model. We follo w closely the dis- 
cussion bv lThompson fc DuncanI (|1993f l. who outline this 
effect for magnetic fields trapped during a core collapse 
of a star leading to a magnetized neutron star. The ba- 
sic idea is that the bubble is threaded by magnetic fields, 
which are dragged along when the bubble expands during 
its rise in the cluster atmosphere. 

Let us consider a magnetic field configuration, in which 
the bubble surface is filled densely with magnetic flux 
loops, analogous to the surface of a star. We further as- 
sume that these loops have a characteristic length scale 
^loopj so that we can approximate their dipole moment 
as /iioop ~ 3$ioop/ioop/87r, where 0ioop is the average 
magnetic flux of each loop. For a random orientation 
of the loops, the net dipole moment of the whole bubble 

can then be written as /if, « -^loop^ioop ~ f'"' ^^Rb,ti.- 
Thus, we can approximate the magnetic dipole field as 
- Bdipoie = ^loopRb.l' a s given by equation (39) of 



[Thompson fc DuncanI ([1993). Assuming that the mag- 
netic field dominates the pressure of the bubble, we can 
then derive the size of a bubble as a function of radius by 
balancing the internal magnetic pressure by the external 
gas pressure. Thus, we can express the bubble radius as 



M 

- Rbfi [1 + [r/vcf] ' . (7) 

with the initial bubble size projected to the cluster center 
given by 



Rb,0 



PO 



(8) 



Note that the exponent in the FMD model is identical 
to that of an adiabatic bubble with an adiabatic index 
of F = 4/3 (model AD43). Thus, one cannot use the 
size evolution of bubbles to distinguish between these two 
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models. Neither can we use this information to constrain 
the fraction of pressure support that is provided from 
randomly oriented magnetic flux loops. 

3.5. Magnetically Dominated Bubbles from 
Current-Dominated Jets ( CD J) 

Three-dimensional magneto-hydrodynamic simula- 
tions of current-dominated jets (CDJ ) have r e cently 

been modeled for the first time (jLi et al.l 120061 : 

iNakamura et all l2006l l2007f) . The jet in these simu- 
lations is launched purely by injecting non-force-free 
poloidal and toroidal magnetic fields. These fields then 
launch a magnetic tower, self-collimate the jet and 
subsonically expand a lobe structure close to the point 
where the background gas pressure starts to quickly 
drop beyond the core radius. 

This jet system is a connected, current-carrying sys- 
tem: the current 1^ first travels along the inner jet axis 
to the lobes, and then returns on the outside of the lobe 
and in a sheath around the jet axis. This results in a mag- 
netic field structure that resembles a tightly wound helix 
inside the jet, with another helical field wrapped around, 
but of opposite sign and less tightly wound. The mag- 
nitude of the 0-component of the magnetic field within 
the lobe is essentially given by applying Ampere's law^: 
= 2c~^Iz R~^- Since the bubbles are magnetically 
dominated, this also sets the size of the bubble which is 
in quasi-pressure equilibrium with the ambient gas pres- 
sure: 

M 

= Rbfi [1 + [r/rc f] ' (9) 

The initial size i?b.o of the bubble projected back to 
the cluster center is then a rather simple function of the 
current Iz fiowing through the system: 

i?M = ^^Po^4 (10) 
V27rc 

To get a better understanding of the magnitude of the 
currents that are expected to be involved, we consider 
a central black hole with a mass around 2 x 1O*M0, 
with a Schwarzschild radius tbh = 2G'AfBH of about 
5.9 X 10"'^'^ cm. Let us further assume an equipartition 
magnetic field of 10'' Gauss in the accretion disk around 
the black hole, that is being threaded through the black 
hole. The current can then be estimated as Iz = 7rr^ j 
with the current density j — (c/47r) V x B. Thus, the 
current^ can be approximated as Iz — (c/4) Br, yielding 
approximately 1.5 x IQ^^ A for a 2 x KrMp, bl ack hole. 
For more details, please refer to lLi et aLl (|2006f ). 

Thus, we can now express equation (jlOp as a function 
of typical values found in galaxy clusters. For a cen- 
tral electron number density of riQ = 10~^ cm~'^ and an 
average temperature of kTo — 5 keV, we get an initial 
bubble size of _Rb.o = 1.021 kpc for a driving current of 
Iz = W^A: 

(11) 

^ We use cgs units throughout this paper. 

^ Note that we used the more commonly used non-cgs unit 
Ampere instead of esus^^ (lA = 2.998 x lO^csu s^^). 



3.6. Summary of Models 

All of the presented models predict that the bubble 
size evolves in a rather simple analytic form: 

Rt=Rb,o [l + ir/rc)T ■ (12) 
Thus, the behavior of the models is very similar, and 
they can be characterized by only two parameters: they 
start with an initial bubble size Rb,o projected back to 
the cluster center, and expand according to the exponent 
a. Table [3] summarizes these parameters for all models. 
Note that Rb^ is related to the energetics of the jet or 
outburst and the central pressure, while a distinguishes 
the nature of the bubbles by the steepness of the evo- 
lution during the rise. Table [3] is organized in increas- 
ing order of a, i.e. the purely hydrodynamic model with 
F = 5/3 (AD53) is the shallowest, while the magnetically 
dominated bubbles in the current-dominated jet (CDJ) 
model expand the fastest. 

4. OBSERVATIONS 
4.1. Bubble Size vs. Radius 

Taking our sample of 64 bubbles in 32 clusters, we can 
now analyze the evolution of the bubble size as a function 
of distance to the center in a statistical sense, to distin- 
guish between the models discussed in the last section. 
Figure [T] shows the bubble size as a function of its phys- 
ical projected distance to the center in kpc. The solid 
line shows the best fit correlation, while the dashed lines 
indicate its la intrinsic width. With an intrinsic width 
of only 0.098, dex, this correlation is surprisingly tight, 
considering the fact that our sample comprises clusters of 
different sizes, richness, black hole masses, etc, and even 
extends into the group and galaxy scale. We also em- 
phasize that this correlation covers almost three orders 
in magnitude and is not simply a consequence of scaling 
both axes with distance. The correlation is as strong and 
as tight when plotted in terms of apparent angular sizes. 
We will discuss in section[5]how incompleteness and other 
effects conspire to produce such a tight correlation. 

As all model predictions scale with [1 + (r/rc)^]", it 
is more instructive to scale the bubble's distance to the 
center by the core radius, as shown in Figure [H The 
triple-dot-dashed line shows the shallowest evolution of 
the AD53 model. The dotted line indicates both the 
FML model, as well as the AD43 model. The steepest 
model predictions shown are the magnetically dominated 
bubbles inflated by current-carrying jets (CDJ, solid line) 
and the continuously inflated hydrodynamic model (CIH, 
dashed line). This diagram is a rather powerful probe, as 
it distinguishes between the power of the outbursts that 
have supposedly created the bubbles, and their subse- 
quent evolution. However, note that the exponent a also 
depends on the /3 parameter, which is centered around 
0.50 for our sample (see Table [1]). All lines drawn in our 
Figures use this average value. 

Figure [3 demonstrates how the AD53, AD43 and FML 
models are all too shallow to represent the data accu- 
rately, while the CDJ model and the CIH model ad- 
equately reproduce the steep slope of the correlation. 
To quantitatively support this impression we produce 
Monte-Carlo simulations of the various models in sec- 
tion [H 

The rather wide spread in the correlation over several 
orders of magnitude in bubble size in Figure [2] must be 
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TABLE 3 

Comparison of bubble size evolution models. Please refer to section [3] for full information on 

THE parameters. 

Model Description ^6 

AD53 Adiabatic expansion of a preformed hydrodynamic bubble with F 
AD43 Adiabatic expansion of a preformed hydrodynamic bubble with F 
CIH Continuous inflation of a hydrodynamic bubble with F = 4/3 

FML Bubble threaded with flux-frozen magnetic loops 
CDJ Current-dominated magneto-hydrodynamic jet model 
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due to the a difference in initial bubble sizes or energy in- 
jection rates. In each of the presented models, the initial 
bubble size and in turn the cavities pV work is related 
to the central pressure pq and the total AGN outburst 
energy, either directly (eq. (O), through means of the 
magnetic dipole (eq. (O) or the current (eq. Thus, 
any spread perpendicular to the indicated lines for any 
model predictions in this figure is primarily due to the 
spread in the initial bubble size, and consequently AGN 
power. 

To test whether the location in the r/vc vs. Rb di- 
agram in Figure [2] is correlated with any other cluster 
characteristic, we have produced color-coded plots with 
different parameters. Two strong correlations emerge, 
which are shown in two top left panels of Figure [H In 
both plots, we have scaled the bubble size by the square 
root of the central pressure, according to equation (|10p 
for the current-dominated bubbles. While this does not 
make a large difference for the general appearance of the 
plot, the y-offsct is now solely due to the current itself 

— 1/3 

for the CDJ model. Scaling the bubble size with p^ or 

Pq as needed by the other models, shows qualitatively 
the same trends. The top left panel of Figure [3] marks the 
data points by central black hole masses, in the sense that 
larger symbols with lighter colors denote more massive 
black holes. Black hole masses are deri ved from veloc- 
ity d ispersion (cr) data from Hv perleda (iPrug niel et al.l 
using the relation given bv lTremaine et al. (200^). 
Thus, this apparent separation of the plot with black 
hole mass could also be rather a separation with a itself. 
The middle left panel of Figure [3] shows an even stronger 
trend with the core radius. Here, the larger symbols and 
lighter colors indicate larger core radii. We will quanti- 
tatively address these trends in section [6] later on. 

In both plots, the data points are mostly separated 
perpendicular to the guidance lines drawn in the plot. 
These guidance lines show a slope of 3/4, the asymptotic 
slope for the CDJ model for a /3 of 0.5, which is the av- 
erage P in our sample. The asymptotic slope of the CIH 
model is 17/24, which is only 1/24 shallower than the 
CDJ prediction. Note that the slope for the next steep- 
est correlation (3/8) is already half as steep as the plotted 
lines, clearly at odds with these models. With a larger 
cavity sample being available in the future, this type of 
study will greatly benefit from splitting the sample into 
appropriate sub-samples that cover only small parameter 
spaces to find out what set of parameters truly governs 
the AGN power. 



As the bubble size have been scaled by the square root 
of the central pressure, the location of the various cavities 
in the plot are solely determined by the current flowing 
through the jet-lobe system - assuming the CDJ model 
is indeed correct. In the middle column of Figure O we 
thus plot the inferred current Iz as a function of black 
hole mass and core radius. One caveat to keep in mind is 
that Iz is derived under the assumption that the geom- 
etry of the observed system is face-on. This projection 
effect will add an uncertainty to the derived current, and 
in general our reprojected bubble sizes and thus our cur- 
rents will be slightly overestimated as a result, on average 
by a factor of 1.6. The projection effect will also add a 
significant amount of scatter into correlations with J^, 
such that 90% of all data points would lie within a fac- 
tor of 2.4. In both plots, the correlation that we find is 
almost linear, as indicated by the dotted guidance lines. 
We find that the inferred currents span about two or- 
ders of magnitude, ranging from ~ 5 x 10^^ Ampere to 
~ 5 X 10^^ Ampere. 

The right panels of Figure [2] show the same analysis for 
the CIH model, but for the energy injection rate instead 
of the currents. The panels show similar trends, with 
only slightly more scatter than for the CDJ model. 

The bottom panel of Figure [2] on the other hand shows 
the same correlations for the integrated 2 cm radio lu- 
mino sity within 1' from the NVSS catalog (| Condon et al.l 
|1998() . The observed trend is much weaker, as can also be 
seen in the middle and right panels. There appears to be 
a large scatter in this correlation, but an apparent lack 
of clusters with high radio luminosity but low current (or 
low energy injection rate) inferred from the cavity sizes. 
This is interesting, as radio luminosity is quite often used 
as a surrogate for total AGN power. This work may sug- 
gest otherwise, at least at the wavelength and the surface 
brightness level of the rather shallow depth of the NVSS 
survey. Radio observations a t longer wavelen gths may 
be better suited for this task (|Birzan 

eEaDHoOSi)- 

4.2. Bubble energy increases with Radius 

It is nowadays standard procedure to measure the size 
of a bubble and the ambient pressure at its location to get 
an estimate of the pV work needed to inflate the bubble, 
which is then in turn used to infer the energy associ- 
ated with the AGN outburst. Interestingly, if one plots 
the inferred energy as a function of radius for the whole 
sample, one finds that the energy increases as a function 
of bubble distance from the center, as shown in Figure |4l 
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Fig. 1. — Bubble size as a function of bubble location in physical 
units (kpc). The light grey data points indicate the data points for 
the multi-cavity system Hydra A. The solid line indicates the best 
fit correlation, while the dashed line show the la intrinsic width of 
0.098, dex. 
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Fig. 2. — Bubble size as a function of bubble location scaled 
by the core radius. The light grey data points indicate the data 
points for the multi-cavity system Hydra A. The lines show the 
predictions from our models for one particular initial bubble size: 
magnetically dominated model (solid lines); purely hydrodynamic 
expansion with an adiabatic index of F = 5/3 (triple-dot-dashed 
line) and F = 4/3 (dotted line); and the flux-frozen magnetic dipole 
model (also dotted line). The grey solid lines show the magnetic 
model for different initial bubble sizes. Note that all lines show pre- 
dictions for deprojected radii, whereas the data show the projected 
radius as measured in the images, as we do not have information 
about the jet geometry. Deprojection will systematically shift the 
data points to the right. Note that the shape and size of the shown 
error ellipse should be also a fair representation of the uncertainty 
in the fitted parameters themselves (red star). 



This trend is obvious in the top panel figure, showing pV 
as a function of the physical radius (in kpc). However, 
there is really no physical reason for bubbles at a larger 
radius to be produced by outbursts of higher energy, the 
radius should simply depend on the rise velocity, and the 
time elapsed since the outburst occured. 

The bottom panel in Figure [3] shows the same data 
points, but this time with the bubble location again 
scaled by the core radius as before. This allows us to 
compare the data to predictions of the various models. 
Since we know the pressure profile as a function of ra- 
dius, and we know the bubble size as a function of ra- 
dius for a given outburst energy, we can easily calculate 
the observable pV work as a function of radius. In fact, 
three out of the five models predict that the inferred pV 
value should decrease with radius. This is due to the fact 
that in the adiabatic model (AD43: dotted line, AD53: 
triple-dot-dashed line) part of the internal energy needed 
to keep up the internal pressure support of the bubble 
is used to expand the bubble adiabatically. In the FML 
moment model (also dotted line), this energy is used up 
to expand the magnetic flux loops on the bubble surface. 
All of these models are again inconsistent with the data, 
along the same lines of argument presented in the last 
section. The only models in which pV is expected to in- 
crease with radius are the CDJ and CIH models. Figure 
[5] shows the correlations between the inferred pV and the 
Iz and Eq values for our sample. Note how a spread of 
only two orders of magnitude in current results in pV 
values spanning almost six orders of magnitude. Part of 
the reason for the intrinsic tightness of these correlations 
is that the parameters are all derived from a common set 
of measured variables (see eq. [10] and [6|) . 

We emphasize that this discussion should not be con- 
sidered an additional argument for the CDJ or CIH mod- 
els, as the pV values are strongly correlated with the 
bubble sizes. We rather want to make the reader aware 
of the severe implications this has on the inferred en- 
ergetics of AGN feedback. It is now accepted standard 
procedure to use the adiabatic model and pV derived 
outburst energies to determine the overall energy bud- 
get of the clusters, and to find out whether AGN heat- 
ing can offset co oling in clusters (e.g. iBirzan et al.ll200 j : 
iBest et aLll2006[ ). However, Figure S] makes absolutely 
clear that we have a very biased view of bubbles in clus- 
ters. If the adiabatic model can really be considered a 
viable model, one needs a way out of this apparent para- 
dox and understand why outbursts seen at larger radii 
tend to have larger pV values. This effect is not only 
apparent in our large sample of different clusters, but 
even within the few individual clusters with known mul- 
tiple bubbles, suc h as Hydra A and Perseus. In Hydra A, 
IWise et aD (|2007[ ) have found a multi-cavity system with 
three bubble pairs at different radii and inferred pV ener- 
gies that continuously increased as a function of radius, 
with the inferred outburst energy se veral times higher 
for th e outermost bubbles. Recently, ISanders fc Fabianl 
(|2007f ) reported the possible detection of another cavity 
in Perseus at very large radii (~ 170 kpc), and also in- 
ferred an outburst energy several times higher than the 
latest outburst. 

There are a few ways out of this dilemma. The first 
possibility is that we simply need a more complex way to 



10 



Diehl, Li, Fryer & Rafferty 




Fig. 3. — The left column shows the scaled bubble radius as a function of bubble location. Colors and symbol sizes indicate the 
dependence on independent third parameters, namely black hole mass (top), core radius (middle) and NVSS radio luminosity (bottom). 
Note how they tend to band the plot along the indicated guide lines, representing the approximate asymptotic slope for the CDJ and CIH 
models. The middle panel plots these same quantities against the inferred current for the CD,1 model. The dashed line indicates a simple 
proportionality between the parameters. Both black hole mass and core radius scale almost linearly with current, though with significant 
scatter. The right panels show the equivalent energy injection rate for the CIH model. Please note that both the current and the energy 
injection rate arc derived from the cavity properties on the assumption of a face-on geometry. Thus, projection issues account for some of 
the scatter in these plots. 
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infer outburst energy from pV , and our current picture is 
too simplistic and thus incorrect. In this case, theoretical 
models have to step up and provide better constraints for 
X-ray observables to see if our simple methods actually 
do have a rigorous basis to stand on. 

The other way out is to assume that nearly all of our 
bubbles undergo continuous inflation, even at this mo- 
ment. However, a large fraction of our bubbles are con- 
sidered "ghost bubbles" and seem to be no longer directly 
attached to the jet or associated with young radio emis- 
sion. On the other hand, many cavities that were earlier 
categorized as "ghosts" now show evidence for older low- 
frequency emission (e.g. lClarke et al.ll2005l ). and m some 
cases an X-ray tunnel of suppressed surface brightness 
may still e xist that con nects to cavities further out (e.g. 
I Clarke et al . 2005; Wi se etal .11200 71 . Creating a purely 
hydrodynamic jet model capable of continuously inflat- 
ing these bubbles at large radii is challenging and - in 
our opinion - unlikely to succeed. In particular, the ge- 
ometric configuration of the cavity system in Hydra A 
(Figure El left panel) would be difficult to explain with a 
conventional hydrodynamic jet model, due to its appar- 
ent "kink" before the outermost cavity pair. A magnetic 
jet on the other hand could be deformed by instabilities 
or ambient gas flow and still be able to deliver energy 
along the field lines to inflate the bubbles even at large 
radii. 

The third and last solution to this paradox is the ef- 
fect of incompleteness, that we are going to address in 
depth in the discussion section of this paper. Under- 
standing this effect will be crucial in correctly addressing 
the question whether AGN can be a general solution to 
the cooling flow problem in clusters and galaxies. To do 
this, we also have to properly understand the size evo- 
lution of the bubbles, as these two effects are of course 
intimately connected. 

4.3. The Multi-Cavity System Hydra A 

Understanding the evolution of bubble sizes in a large 
sample of clusters gets complicated by the fact that the 
sample is a combination of a wide variety of cluster prop- 
erties. A much cleaner approach is to look at the evo- 
lution of bubbles in a single system. Since we are of 
course unable to observe a single cavity multiple times 
during its rise in the cluster due to the extremely long 
rise times (several Myr) , we have to rely on systems that 
harbor multiple cavities at different distances. Based on 
the reasonable assumption that AGN produce at least 
similar, if not identical bubbles over time, we can then 
use these individual bubbles to infer the size evolution of 
a rising bubble over a ran ge of distanc e s. 

In the case of Hydra A, IWise et afl (|2007f ) lately dis- 
covered a giant multi-cavity systems consisting of 6 dis- 
tinct cavities, with the outermost bubble being intact 
at a distance of over 200 kpc from the center. The left 
panel of Figure P show s the cavity system, reproduced 
from IWise et al.r(|2007l ) with permission from the au- 
thors. The properties of these bubbles are also listed 
in Table H 

Since Hydra A is now the best multi-cavity system with 
adequate published data to test our models, we have also 
marked its data points with light grey in Figures [l] [2] 
and m The right panel of Figure [5] shows the cavity sizes 
as a function of distance from the center for only the 
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Fig. 4. — Top panel: Bubble pV energy as a function of physical 
radius. Note that there is a trend of bubbles in the outskirts of 
clusters to have higher energies. The light grey shaded data points 
show the multi-cavity system Hydra A. Bottom panel: Same plot 
but with the radius scaled by the core radius of the cluster. The 
lines show the predictions from our models: the CD J (solid line), 
CIH (dashed line), AD53 (triple-dot-dashed line) and AD43 (dot- 
ted line), and the FML model (also dotted line). Note that only the 
CDJ and CIH models actually predict an inferred increase in pV 
energy with radius, all other models predict a declining behavior. 



Hydra A cavities. Note how quickly the cavities expand 
as they rise outward. They tend to grow much faster than 
expected for the purely hydrodynamic models AD53 and 
AD43 (triple-dot-dashed and dotted lines, respectively) 
and the magnetic dipole model (FML, also dotted line). 
The only two models in agreement with all six cavity sizes 
are the CDJ and CIH models again. If the CDJ model is 
indeed correct, the jet in Hydra A carries a current of ^ 
1.5 — 2 X 10^^ A, for the CIH model, this would require an 
average energy injection rate of '--^ 0.8 — 3 x 10'*'* erg sec~^ 
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Fig. 5. — Bubble pV energy as a function of inferred energy in- 
jection rate for the CIH model (top panel) and current for the CDJ 
model (bottom panel). Note that in each plot, both parameters are 
derived from common parameters, which may account for some of 
the tightness in the correlation. However, the correlation with the 
current is significantly tighter than with energy injection rate. The 
remaining scatter is consistent with projection effects alone. 

over the largest cavity life time of ^ 230 Myrs. 

A look at the radio emission of Hydra A shows that 
these cavities all seem to be connected at lo w radio fre- 
que ncies out to a few hundred kiloparsec (jLane et al.l 
|200 4). There also seems to be "kink" in the radio emis- 
sion close to the location of cavities C and D (Figure [6]). 
This may suggest that the outer cavity pair has a differ- 
ent orientation than the other two pairs, and that it has 
a larger inclination angle with respect to the plane of sky. 
This is consistent with the data points for C and D in 
the right panel of Figure [5] to be shifted towards smaller 
radii due to projection effects. This would also explain 
why the inferred current (or energy injection rate) for 



the outermost bubble pair is slightly larger than for the 
two inner bubble sets. 

For projection effects to cause individual disconnected, 
purely hydrodynamic bubbles to produce this steep plot, 
the jet would have to conspire and be significantly closer 
to the line of sight with each successive cavity pair. While 
this is sufficiently unlikely in itself, the outermost bubbles 
would also have to be several times further away when 
deprojected, at which point the contrast in the bubble 
would be way too small to be detectable due to the drop 
in surface brightness at larger radii. 

4.4. The Multi-Cavity System Perseus 

The other well-known mult i-cavity systems is the fa- 
mous Perseus clusters (e.g. iSanders k, Fabiaiil [2007.) . 
Judging from t he appearance of t he inner set of bub- 
bles in Perseus ([Fabian et al.ll200^ . the geometry of the 
cavity set right next to the cluster center seems to be 
rather complicated, and different from that of the outer 
set, both of whi ch are part of our current data set. 

Only recently, ISanders fc Fabiaiil ([2007[ ) noted a possi- 
ble intact cavity on the outskirts of the Perseus cluster 
around 170 kpc. This tentative cavity candidate has a 
width of approximately 73 kpc (200", Jeremy Sanders, 
private communication). Unfortunately, the bubble is 
located off the main ACIS chip, and is thus not included 
in all exposures. In addition, the Chandra point spread 
function is rather large at this location, and the surface 
brightness level very low, which makes its detection diffi- 
cult. As a result, the height of the possible cavity remains 
very uncertain, and can be anywhere from around 9 kpc 
(25") to 25 kpc (65"), depending on the surface bright- 
n ess model that one choose s to subtract. 

ISanders k. FabiaiJ ([2007D already note that the out- 
burst energy associated with this new bubble requires an 
outburst of at least twice of that needed to create the 
inner bubbles. As we have pointed out in the last sec- 
tion, this could also be possibly explained by the CDJ or 
CIH models. Figure [7| shows the size evolution of the 4 
observed cavities in the Perseus cluster (black filled cir- 
cle). For the newly detected cavity candidate, we show 
two estimates for its bubble size: The red circle shows 
the width as the bubble size, the green circle uses the 
geometric average between the width and the smallest 
height estimate. The lines flowing through these two 
data points indicate how the bubble would have evolved 
to its current size assuming our different models. Unfor- 
tunately, our results remains inconclusive. If we assume 
that the green data point is representative of the bub- 
ble size, the energetics of the outburst that created this 
bubble is comparable to the outburst that is currently in- 
flating the innermost set of cavities. However, if we take 
the red data point as given, the adiabatic models yield an 
outburst roughly double as energetic. It is intriguing to 
note that both the CIH and CDJ models are consistent 
with the larger bubble size estimate. A deep observa- 
tion pointed directly toward this outer bubble (to have 
a smaller Chandra point spread function) is most likely 
needed in order to firmly assess the true size of this bub- 
ble and to paint a better picture of the outburst history 
of Perseus. More deep observations of other multi-cavity 
systems are needed to further distinguish between these 
models. 
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Fig. 6. — Left: The multi-cavity system in Hydra A, reproduced from 'Wise et alJ (2007') with permission from the authors. The black 
area is excess X-ray emission left-over after an elliptical surfa ce brightnes s mod el has been subtracted. Right: Data Points: Bubble sizes for 
Hydra A as a function of distance to the center, taken from lWise et al] II2007I ): Lines show predictions from the AD53 (triple-dot dashed 
line), AD43 (dotted line), FML (also dotted line), CIH (dashed line), as well as the CD J model (solid line). The cavity labels are the same 
in both plots. 
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Fig. 7. — Bubble sizes for Perseus as a function of distance to the 
center. Lines as in Figure |6] The red data point shows the upper 
limit for the new bubble size estimate, the green data shows a lower 
limit. The correct answer will likely lie somewhere in between these 
two extremes. 

5. DISCUSSION 

5.1. Projection Effects 

Projection has several effects on the observables that 
are crucial to our analysis. The most obvious effect is 
that we are measuring the projected distance to the cen- 
ter, rather than the true distance. This will systemat- 
ically shift the data points toward smaller radii in all 
plots. In fact, all of our radii should be considered lower 



limits to the true location of the bubbles. This will not 
only affect the radii themselves, but also the point at 
which other quantities are evaluated at, like density, tem- 
perature and pressure. In general the temperature rises 
outward in these systems, thus the temperature at the 
location of the bubble is likely to be systematically un- 
derestimated. The density and ambient pressure on the 
other hand will always be overestimated. This also means 
that any rise times derived from using the projected ra- 
dius rather than the true distance to the center will result 
in estimates for the rise times that are systematically too 
low. We also note that the smaller the observed radius 
is, the higher the probability that it is due to an effect 
caused by projection. 

But there are more subtle effects that projection has on 
our data. As we do not have an automated tool to detect 
bubbles, one has to rely on human experience in finding 
and identifying these systems. This task is much more 
difficult, if the cavities overlap with the bright cluster 
center or the bubble on the opposite side of the cluster. 
In fact, our sample does not contain any cavity system in 
which the bubble size exceeds the projected distance to 
the center, the slope of which is shown by the black solid 
line in Figure [51 even though this is statistically very 
improbable. This suggests that our sample is affected by 
what we will refer to as a "geometric" selection effect, 
introduced by our manual detection process. 

5.2. Incompleteness Effects 

One very important issue in a systematic study of 
X-ray cavities is to address incompleteness, i.e. how 
likely it is to detect a bubble of certain properties. 
lEnfiiin k Heind ()2002t ) address this issue by computing 
analytical formulae for the S/N ratio of a spherical hole 
in the isothermal /? model. In this context, the signal 
5* is the deficit in counts due to the hole, i.e. the sum 
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Fig. 8. — Grey shaded data points are the same as in Figure 
[2] The colored Unes shows the required size evolution of a butble 
to keep a const ant signal to noise ratio while rising, taken from 
IEni3hn fc H"em3 ll200a') . for different angles of inclination with the 
plane of sky; 0° (face-on, red), 30° (green), 60° (dark blue), and 
80° (light blue). The dotted red lines repeat the line for face-on 
detectability for different initial bubble sizes, spaced by a factor of 
4. The black solid line indicates a slope of "1" , to show equality be- 
tween bubble size and projected radius, our "geometric" detection 
criterion. 



of counts that would be emitted in the cavity region ac- 
cording to the emissivity profile. The noise N on the 
other hand is defined as the square root of the sum of 
counts from the surrounding regions contributing to the 
projected area that the cavity occupies. Enfilin & Heinz 
find that the S/N of an observed bubble depends strongly 
on their location in the cluster, as well as the specifics of 
the observations. Using our nomenclature, we generalize 
their equation (11) for any bubble of size at a radius 
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(13) 

where 6 is the angle by which the jet axis is tilted out of 
the observed sky plane, 5 is the ratio between the peak 
flux and the background flux (typically 5 «1, but this 
depends on observational details), At is a function of the 
detector's effective area and the distance of the target, 
and B is the incompl ete beta function. Please refer to 
lEnBlin fc HehI3 ()2002[ ) for further details. 

Figure [8] shows how the size of a bubble needs to evolve 
in order to keep a constant S/N throughout its lifetime, 
i.e. to be equally likely to be detected. The colored solid 
lines show the same line, but for different angles with 
respect to the line of sight. The dotted red lines repeat 
the line of constant detectability for the face-on case for 
different initial bubble sizes i?b,0; spaced by factors of 
4. This plot dramatically emphasizes the importance of 



understanding incompleteness for our sample, as the data 
points (grey shaded) mostly follow a similar slope as the 
incompleteness limits. We would like to stress this under- 
appreciated fact and urge other researchers to take this 
into account when drawing conclusions about the energy 
budget of AGN in clusters. 

5.3. Hydrodynamic Instability Effects 

One of the most astonishing discoveries about cavities 
in clusters of galaxies has been that these bubbles seem 
to stay intact over a very long time. There are a num- 
ber of hydrodynamical instabilities that should develop 
in these rising bubbles, namely Kelvin-Helmholtz (KH), 
Rayleigh- Taylor (RT) and Richtmyer-Meshkov (RM) in- 
stabilities. Here we use simple estimates of the growth- 
time of these instabilities to both determine whether such 
instabilities will produce observable features in the cav- 
ities and determine which instabilities are most impor- 
tant. 

As the bubble rises, the shear flow between the 
low-density rising bubble and the denser surrounding 
medium will develop Kelvin-Helmholtz (KH) instabili- 
ties. A simple estimate of the growth time for such an 
instability can be obt ained from perturbation theory (e.g. 
lChandrasekhMl[l96ll) : 
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where A/i is the size scale of the instability and Au is 
the velocity of the rising bubble relative to its surround- 
ing medium. This bubble velocity can be approximated 
by using several estimates, such as the terminal buoy- 
ancy speed, the r efill time, or the s ound speed in the 
ambient medium (|Birzan "etan[200l . In d we decide 
to list the inte g rated buoyant rise times according to 
lEnBlin fc Heind ()2002l ) for an adiabatic expansion of the 
bubbles, pb and pamb are the densities of the rising bub- 
ble and its surroundings respectively. Typically, values 
of pb < 0.1 — O.Olpamb are to be expected from numerical 
simulation and from limits on the X-ray emission coming 
from these cavities. In our estimates, we assume a fixed 
contrast of 0.01 between the densities, and a rise veloc- 
ity close to the sound speed. In Table O we list the KH 
igrowth times for the largest mode that could possibly 
develop, i.e. for Xh = Rh, the current size of the bub- 
ble. As the density contrast is rather large, the growth of 
KH instabilities should be rather slow compared to the 
inferred rise times of the bubbles. Assuming a density 
contrast of only 0.1 instead decreases the growth times 
by a factor of ~ 8, which would bring them closer to the 
vicinity of the bubble ages. Thus, we expect KH insta- 
bilities to only play a role on a scale much smaller than 
the current bubble size, at least A/i < 1/10 i?b- While it 
is unclear what the exact observational signature in X- 
rays should be, Chandra^s resolution should in principle 
already be good enough to find signatures of developing 
KH instabilities. The fact that the sides of the bubbles 
seem to be completely intact calls for either a suppres- 
sion of KH by viscosity of magnetic fields, or a very high 
density contrast. 

The top of the bubble on the other hand is prone to de- 
velop Rayleigh- Taylor instabilities, as a heavy fluid (the 
ambient medium) is accelerated on top of a light fluid 
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(the bubble), with the gravitational acceleration point- 
ing toward the light medium (toward the center of the 
cluster). Estimates of the growth time for RT instabili- 
ties are often given as: 



RT 
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where g is the gravitational acceleration and A is the at- 
wood number: (pamb — Ph)/ (pb + Pamb) • Because the den- 
sity contrast is supposed to be high (pb < 0.1— O.Olpamb), 
the atwood number'' is not larger than about 1.2. Ta- 
ble O lists our RT growth time estimate for an atwood 
number of 1.0 and the largest possible mode, i.e. the 
current size of the bubble. The gravitational acceler- 
ation has been estimated using the /3 model assuming 
hydrostatic equilibrium. While these growth time esti- 
mates are thus not intended to be used at face value, they 
show one characteristic that is ubiquitous for the whole 
sample: they are all on the order of or smaller than the 
rise time of the bubbles. This is particularly interesting, 
since we listed the growth time for the largest possible 
scale length, which is the longest. Thus, RT instabili- 
ties should at least have already begun to shred most of 
the bubbles apart on scales of a significant fraction of 
the bubble size. In fact, the RT instabilities would have 
had sufficient time to grow a global instability on the 
bubble siz e scale for most o f the bubbles (Figure [H blue 
squares). iPizzolato et al.l (|2006f ) have argued that the 
one should take into account the expansion of the bubble 
during inflation when computing the RT time scale. This 
acceleration diminishes the effect of the gravitational ac- 
celeration, leading to a suppressed onset of RT instabil- 
ities. However, this theoretical work is rather idealized 
and it has yet to be determined if this mechanism actu- 
ally works successfully in realistic 3-dimensional hydro- 
dynamic simulations. We have only considered the tra- 
ditional Rayleigh- Taylor time scales but keep this caveat 
in mind. 

However, probably the fastest growing instability for 
the case of our rising bubbles is the Richtmyer-Meshkov 
(RM) instability. RM instabilities grow when low- 
density material (the cavity ) "plows" through a higher- 
densi ty ambient medium (jRichtmveil 119601 : iMeshkovl 
Il969| ). These instabilities are usually important when 
the contact surface between two fluids is accelerated (e.g. 
by a shock passing through). The characteristics of this 
instability are usually the generation of "spikes" sepa- 
rated by voids on the surface. An analytic estimate of 
the RM instability is given by: 



(16) 



This instability is the fastest growing (at all size-scales) 
of the three discussed in this section, as shown in Tabled 
Again, here we used the bubble size as the largest possi- 
ble mode, and the sound speed as a proxy for the cavity 
velocity. Note that the RM growth times are all on the 
order of Myrs or below, always either on the order of the 
bubble rise times, or up to an order of magnitude shorter 

Note that the more general estimate for the RT growth time 
requires knowledge of the entropy gradient at the instability sur- 
face. But for our estimates, the simple atwood number prescription 
provides a safe upper limit on the growth time. 
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Fig. 9. — Adiabatic rise times for an assumed F = 4/3 against 
the instability time scales for the largest modes (A^ = ij^) of the 
Richtmyer-Meshkov (RM, red circles) and Rayleigh- Taylor instabil- 
ities (RT, blue squares). The solid line delineates equality between 
the two timescales. Every bubble above this line should have al- 
ready developed instabilities on the order of the bubble size. The 
three dotted lines (top to bottom) show equality smaller modes, 
namely \h/Rh = 1/2 and l/^/8 for the RT instability and 

Ah/_Rb = 1/2, 1/4 and 1/8 for the RM instability. 



than them (Figure[51 red circles). Thus, we conclude that 
RM instabilities almost certainly should have developed 
on at least part of the top surface of the bubble. 

We urge the reader to remember that these timescales 
are approximate. It may be that although we predict 
that instabilities on sizescales equal to the bubble size 
should develop for most of our cavities, uncertainties in 
such growth-time estimates prevent us from stating this 
with complete certainty. However, the predicted growth 
time for RM and RT instabilities for 1/lOth the bubble 
size is always over a factor of 10 lower than the rise time 
and such instabilities almost certainly should have devel- 
oped (Figure [9]). Such instabilities should be observable 
with the high resolution of Chandra and are not seen to 
our k nowledge, ev en in the very deep observations of Hy - 
dra A I Wise et all (|2007|) or Perseus (|Fabian et al.ll2006f) . 

RT and RM instabilities should be entirely suppressed 
for magnetically dominated models for modes larger than 
the coherence length of the magnetic field, which for the 
CDJ model is on the order of the bubble size. 

Several caveats should be kept in mind when using 
these simple estimates. First of all, the rise times that we 
have computed assume that the bubble is inflated at the 
center. If this assumption is far off, and the bubbles get 
inflated at radii on the order of the core radius, then the 
quoted rise time should be considered an upper limit on 
the rise time. If this is the case however, the jets should 
be highly supersonic, which would result in the shocks 
and hot spots at the cavity rims in general. Projection 
effects on the other hand will cause the rise times to be 
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underestimated. Which of the two effects dominates or 
whether they cancel each other out on average, is unclear 
and depends on the details of the jet dynamics. 

6. MONTE CARLO SIMULATIONS INCLUDING 
INCOMPLETENESS AND PROJECTION 

6.1. Methodology 

To properly address the question of which model can 
truly reproduce the properties of the cavity sample, one 
would have to ideally model the complete manual detec- 
tion process, a very work intensive procedure. We follow 
a somewhat simpler approach by using a Monte-Carlo 
(MC) technique that mimics the detection process with- 
out the need of human intervention. The outline of the 
MC process is as follows: 

(i) Pick the initial bubble size Rb^ (i.e. the outburst 
energy) 

(ii) Pick the distance r to the center and pick (p and 6 
to randomly orient the bubble with respect to the 
line of sight 

(iii) Determine the probability of detecting the cavity 

(iv) Repeat steps (i)-(iii) until you have "detected" as 
many cavities as in the original cavity sample (64) 

(v) Fit the simulated r-Ri, data set to determine the 
parameters slope, y-offset and intrinsic width 

(vi) Repeat steps (i)-(iv) 1000 times to find the uncer- 
tainties in the parameters 

The key challenges in this Monte-Carlo bootstrap 
method are in picking the correct initial bubble size (step 
i) and in determining whether the evolved bubble would 
be detected or not (step iii), which are unfortunately also 
the steps that are the most uncertain. During the rest of 
this section, we we will follow this basic methodology, but 
vary the details of these steps. For a valid model, the MC 
simulations should be able to reproduce the three main 
characteristics of the distribution of data points in the 
r-Ri, plots, namely the slope, y-offset and the intrinsic 
width, as de termined by t he fitting method bandf it (see 
appendix in lDiehl fc Statlerl (j2008af l for a description of 
this method). 

The amount of work that has been done on the de- 
tectability of X-ray cavities is very limited. We employ 
the following criteria to determine whether a bubble is 
detected or not: 

(a) S/N criterio n: We co mpute the S/N of a bub- 
ble from Enfi lin fc Hein z (2002) using our equation 
(fT3|) . For simplicity, we then use a Gaussian prob- 
ability distribution to compute the detection prob- 
ability. 

(b) Geometric criterion: If the extent of the bubble 
overlaps with the center of the cluster, the cavity 
would be extremely difficult to detect, in particu- 
lar if there is another bubble on the opposite side 
of the cluster that overlaps then as well. In prac- 
tice, we cut every bubble whose radius exceeds the 
projected distance to the center. 

(c) Hydrodynamic Instability: In the purely hydrody- 
namic models, we remove all cavities whose rise 
times exceed the RT or RM growth times accord- 
ing to equations ([15]) and p6|) . However, tests show 



that this has very little to no influence on the sim- 
ulations. 

(d) Detector Size: If the observed distance to the cen- 
ter exceeds half the size of the Chandra CCD (4'), 
we remove the candidate from the sample, as most 
observations in our sample are single-CCD obser- 
vations. 

For more details on any of these criteria, please refer back 
to section [5] 

Thus, this process is rather complex and depends a 
large variety of input parameters. One key parameter 
is the assumed shape of the underlying pressure profile, 
given by the /3-profile parameters P and Tc and po (i.e. 
no and kTg). The steepness of the pressure gradient is 
crucial in determining how quickly the bubbles expand 
while rising. 

Further, observation specific details become important 
for the detection probability (see section [521 for more de- 
tails): the effective exposure time r, the average detector 
efficiency e, the X-ray background level Shg, as well as 
the distance to the object d. In the following discussion, 
we will describe our choices for all of these parameters 
and lay out how these choices influences the behavior of 
the simulations. 

6.2. Basic Simulations 

We start with a MC simulation with the most simplis- 
tic assumptions possible and demonstrate why this is in- 
sufficient to reproduce the observations. The most basic 
approach one can take to design a replica of a true sample 
to O*'' order is to use one set of fiducial template cluster 
properties common to all 64 cavities. In these MC simu- 
lations, the morphological and observational parameters 
of this one template cluster are chosen as the average of 
the cluster sample, and kept fixed for all cavities. 

It turns out that no model comes even close in repro- 
ducing the cavity sample's properties. Slopes are gen- 
erally too shallow. In fact, the simulated samples be- 
have mostly as expected from the theoretical considera- 
tions. The hydrodynamic models AD53 and AD43 and 
the FML model show the shallowest evolution, but get 
steepened due to the incompleteness cutoff for small cav- 
ities at large radii. The CD J and CIH models are affected 
the least by this incompleteness effect, as the evolution 
is essentially parallel to the incompleteness cutoff, i.e. a 
bubble that would be detected at small radii is likely to 
be detected further out as well. 

Since none of the simulations can reproduce the obser- 
vations, we have to add additional sophistication to our 
MC simulations. As our sample encompasses a very large 
variety of environments, pressure profiles, size scales, 
temperatures, etc., we need to take the intricate inter- 
dependence between these parameters into account (e.g. 
P with kTo). As we will show in the following sections, 
the complex correlations between these parameters are 
the reason this simple test fails. However, we empha- 
size that the differences in the predictions of the various 
models are the strongest in this simplified case. Thus, 
this type of analysis will greatly benefit from having a 
larger cavity sample in the future, that can be split up 
into groups with very similar properties. We tested this 
idea with the current sample, but its size of 64 objects 
does not allow for sufficient statistical accuracy on the 
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sub-samples to constrain any models. 

6.3. Use of a Full Set of Parameters 

In order to improve on the basic simulations, we have 
to use as much of the cluster properties as possible to re- 
produce realistically the sample properties. Since we do 
know the details of the cluster profile and the observa- 
tions, we use all of the available parameters to simulate 
our set of cavities: (3, rc, kTp, tiq, pq, d, r, e, and S'bg. 
Thus, for each cluster in our sample, we can choose an 
initial bubble size and simulate cavities specifically tai- 
lored to each cluster and its observation as well. 

The observed sample itself is essentially a random con- 
glomerate of observations of clusters, groups and galaxies 
for which cavities have been detected. To recreate this 
complex selection function, we decide to simulate cavities 
individually for each object and repeat until the exact 
number of observed cavities has been matched. 

The advantage of this approach is that it is very tightly 
constrained, and that we take our host selection function 
properly into account. The only free parameter in our 
MC simulations is now the choice of our initial bubble 
size, and we explore several different prescriptions in how 
to choose this parameter. First we use the measured cav- 
ity sizes to infer back the initial bubble size based on the 
cavity model that we are investigating at the moment, 
i.e. we use our theoretical understanding of bubble size 
evolution (equations (O, (O, and ([9|)) to calculate the 
size the cavity would have started with at the center. 

The only problem with this method is that we have 
to assume a geometry, as we only have the projected dis- 
tance to the center available, and that the observational 
errors in determining the observed bubbles will affect the 
bubble sizes in our simulations. It turns out that this par- 
ticular problem makes it impossible to reproduce the full 
observed sample properties. Figure fTOl shows the results 
of the linear fit to our simulated equivalent of Figure 
[1] {Rb vs. r). The plot shows the 3d-parameter space 
(slope, y-offset and intrinsic width) occupied by our sim- 
ulated models (light blue: AD53, dark blue: AD43 and 
FML, orange: CIH, green: CD J), along with the ob- 
served parameters. As one can see, none of the models 
agree with the data (red star). The current-dominated 
model comes closest, with the parameter set being just 
outside the 3(T contour lines. The main problem with the 
models in general is that the extremely tight width of the 
correlation cannot be reproduced. We think that this is 
the result of the amplifications of the observational er- 
rors and the projection effects in the actual sample. In 
fact the difference between the inferred intrinsic widths 
of the simulated samples to the observations is close to 
the value expected from observational error alone (a 10% 
relative error is 0.4dex). The error in the initial bubble 
size will not only affect the intrinsic width, but also slope 
and offset of the correlations. 

The fact that the observed correlation between r and 
Rb (Figure [1]) is so tight, suggests that the initial bubble 
size Rb.o is intimately connected to other intrinsic cluster 
properties, which we explore in the next section. 

6.4. Initial Bubble Size Rb^ Linked to Core Radius rc 

As the outcome of the last set of MC simulation using 
the full set of sample parameters has shown promising 
improvement over the basic simulation, we now look for 



the best way to set up the initial bubble sizes. Figure set 
[3] suggests two possible parameters: the black hole mass 
TTiBH, and the core radius rc- 

Using the black hole mass as a surrogate for Rb^o does 
not substantially improve the results of the last MC 
simulation, though it does move the simulation closer 
to agreement. However, a substantial improvement is 
achieved when relating Rb^ to the core radius. We ex- 
plore various power-law correlation, but find that a sim- 
ple proportionality of the type Rb^o = ^rc works best. 
Figure [TT] shows the parameter space for this MC simua- 
tion. Note how only the CDJ model is able to reproduce 
the data, while all other models can be excluded at more 
than the 5a confidence level. 

We would like to emphasize that even with a large 
number of combinations of the parameters available to 
us, we did not manage to bring either the hydrodynamic 
or the magnetic dipole model into agreement with the 
observed data, even when only using subsets of the data. 
In particular, even with fine-tuning, not a single MC sim- 
ulation for the AD43, AD53, FML or even CIH models 
has been able to produce a correlation as tight as the 
observed correlation. 

Figure [T^] shows a sample of 5 of our MC simulations 
for the various models (left to right: CDJ, CIH, FML 
and AD43, AD53). These samples have been drawn ran- 
domly, and not been hand-picked by us. Note that the 
CDJ samples closely resemble the observed sample (bot- 
tom, for reference), followed by the CIH model, while the 
other two columns do not. They exhibit signs of "band- 
ing" in the center, where one can actually see the effect 
of the shallower slope for the evolution. This effect not 
only makes the fitted slope shallower, but also widens the 
correlation significantly, which in turn affects the inferred 
intrinsic width and the Y-offset. 

We conclude that with our current understanding of 
incompleteness issues, the CDJ model is the only model 
that we considered that can quantitatively explain the 
cavity data. Any other correct cavity or jet model would 
have to mimic the behavior of the CDJ model in terms 
of the size evolution of cavities. We also explored the 
possibility of an unknown model whose cavities expand 
even faster than that of the CDJ model, by simply in- 
creasing the exponent a. It is interesting to note that 
those models also fail to reproduce the data. The CDJ 
model seems to lie in the "sweet spot" for the observed 
tight r — Rb correlation to work. 

Surprisingly, even the CIH model can be formally ex- 
cluded by the MC simulations, even though its predic- 
tions are close to the CDJ model at large radii. The 
differences at small radii or the slightly shallower slope 
at large radii (~ 17/24 instead of ~ 3/4) are sufficiently 
to be a worse match. This trend is supported by the re- 
lations of pV with Iz and Eq in Figure EJ where the CDJ 
model exhibits a much tighter correlation. Nevertheless, 
the predictions for the CIH model are sufficiently close 
to the CDJ model that details in the cavity detection 
algorithm may yet change this answer. 

6.5. Caveats 

While the above simulations provide good constraints 
on the nature of cavities in clusters, we have to be care- 
ful in not over-interpreting their results. The MC simu- 
lations are based on the assumption that we in fact do 
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Fig. 10. — Three-parameter plot for one of our Monte-Carlo simulations of section [6.31 We produced 500 simulations for each model 
with 64 detected cavities each. The plot shows the best linear fit parameters: slopes, y-offsets and intrinsic widths. In this run, the initial 
bubble size was derived from the observed bubble sizes. Observational errors and projection effects results in an additional scatter that 
increases the width of the observed correlation, preventing the simulations to match the observations. Colors indicate the model: CDJ 
model (green circles); CIH (orange diamonds); FML (dark blue triangles); AD43 (also dark blue triangles) and AD53 (light blue boxes). 
The red star indicates the location of the best fit through the actual data, that we want to reproduce. The contours show the 1, 2 and 3(t 
confidence levels for the CDJ model. Note that the CDJ model gets closest to the observed correlation, though still being slightly more 
than 3cr away. The shape and size of the shown error ellipse should be also a fair representation of the uncertainty in the fitted parameters 
themselves (red star). 



understand all the incompleteness issues associated with 
the data. However, this is not really the case. While 
lEnBlin fcT leinz (2002) provide a nice analytic estimate 
of the SjN of a cavity in a cluster atmosphere, it is un- 
clear how this estimate translates into the real world of 
actually identifying the cavities. This task is left to expe- 
rienced astronomers looking at the data individually and 
detecting bubbles by eye. The only tools guiding them at 
this task are sometimes radio observations and unsharp- 
masking. The presence of full or at least partial cool 
rims around cavities greatly enhances their detectability, 
which may bias our view of cavities in general. 

There is no automated tool available to detect bubbles. 
Until this is available, and we can conduct quantitative 
tests of our ability to detect bubbles, we cannot claim 
that we truly understand the incompleteness effect for 
cavities in clusters. 

The fact alone that we have to assume an ad-hoc ge- 
ometric detection criterion shows how fragile this con- 
struct is. Our MC simulations are very sensitive to in- 
completeness issues, which are essential in producing the 
tight correlations seen in the data. We suggest that the 
type of analysis described here should be repeated when 
incompleteness is better understood and a larger sample 
is available. 



7. CONCLUSIONS AND FUTURE WORK 

We present the first systematic analysis of cavity sizes 
from cavities located in cluster, group and galaxy envi- 
ronments. We summarize predictions from four partic- 
ular models: three purely hydrodynamic (AD53, AD43 
and CIH) and two magnetically dominated models (FML 
and CDJ, see Table [S]). The two simplest hydrodynamic 
models AD43 and AD53 assume the bubble is formed a 
priori close to the center and its interior pressure due to 
gas pressure from a hot, low density gas. The bubbles 
then adiabatically expand during the buoyant rise in the 
intracluster medium. These two models differ only in the 
use of two different adiabatic indeces F = 5/3 and 4/3. 
The CIH model also assumes adiabatic expansion, but 
rather continuously injects energy into the cavities. This 
requires the jet to stay connected to the cavities as they 
rise, which may be challenging to achieve in reality. 

The first magnetically dominated model (FML) as- 
sumes that the magnetic field configuration consists of 
a dipole magnetic field generated by randomly oriented 
magnetic fiux loops. We show that the predictions from 
this model are identical to the F = 4/3 hydrodynamic 
model, which makes a distinction of these models from 
cavity sizes alone unfeasible. Thus, one cannot possi- 
bly use cavity sizes to determine the fraction of pressure 
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Fig. 11. — Same plot and data symbols as in Figure [TOl but this time for the simulation described in section [6.41 In this run, the initial 
bubble size is set to be proportional to the core radius. The contours show the 1, 2 and 3a confidence levels for the current-dominated 
model, with the la contour already encompassing the data parameters in all three plots. Note that all other models can be excluded on 
more than the 5a level. 



supplied by magnetic loops threaded through the bubble 
surface. The last model that we consider is a magneti- 
cally dominated model in which the jet carries a current 
(CD J). In this model, the bubble size is set by the point 
where the interior magnetic pressure due to the flowing 
current is matched by the outside pressure. 

Our analysis of the 64-cavity sample shows that the 
bubble sizes are much bigger than predicted from the 
adiabatic and the magnetic dipole model, and that they 
are more consistent with the continuously inflated hydro- 
dynamic model and the current-dominated jet models in- 
stead. In particular we find that the inferred currents in 
these systems scale almost linearly with the core radius 
of the cluster and the central black hole mass. 

We conduct Monte-Carlo simulations of the bubble de- 
tection process including incompleteness issues, which we 
also discuss in length. We show that the magnetically 
dominated model is always favored compared to the other 
four models. For the case that the initial bubble size is 
linked linearly to the core radius, the current-dominated 
model is consistent with the data within one standard 
deviation. Even with fine-tuning, we are unable to bring 
the other four models into quantitative agreement with 
observations, and never get closer than a, 5cr offset. 
Any viable future jet model has to at least mimic the 
behavior of the CD J model in terms of the size evolution 
of bubbles. 

In order for the purely hydrodynamic model to work, 
we also need a mechanism to efficientl y suppress instabil- 
ities. Viscosity (|Revnolds et al.|[2005D or magnetic drap- 



ing (jRuszkowski et al.]|2007D may provide this added sta- 
bility. Continuous inflation may suppress the growth 
of Ray leigh- Taylor in stabilities on top of the bubble 
(jPizzolato et al.|[2006f ). A fine-tuned jet model continu- 
ously inflating cavities even at large radii may also result 
in a faster size evolution of the cavities. However, bring- 
ing such a model into agreement with constraints from 
both X-ray observations and old particle ages inside the 
cavities will be challenging. 

We stress the fact that we do not consider our results 
as final. One of our key assumptions is that the bub- 
bles are always in pressure equilibrium with the ambient 
gas. We base this assumption on the general lack of of 
shocked gas surrounding the bubbles. If this is incorrect 
and the bubbles are overpressured closer to the center, 
they would expand faster than the simple AD53 or AD43 
models that we have considered here. To be a bit more 
quantitative, let us take a fiducial observed bubble of 
size i?;, at r = 2 Tc in a cluster with = 0.5 and as- 
sume that it is currently in pressure equilibrium. If we 
further take the steeper slope of the CDJ evolution as 
a given, we can reproject the bubble size back to the 
center: Rbfi = i?fc [1 -f (r/rc)^]"^/*^ « 0.55 Rt- Since we 
know its pressure and volume now, as well as its size be- 
fore, we determine the the bubble's internal pressure pb^ 
it must have had to start with: pb^o = Pb (Vb/Vb.o)*^^ = 
Po [1 + (r/re)2]-3/4 {Rb/Rb,o)\ assuming T = 4/3. This 
yields an overpressure oipb^/po ~ 3.34. Such an over- 
pressure would have driven a shock that should have been 
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Fig. 12. — Sample plots from our Monte-Carlo runs. The bottom panel shows the original data from our cavity sample. The columns 
shows 5 of our Monte-Carlo runs for the CDJ, CIH, AD43/FML, and AD53 models (from left to right). The green lines in each plot show 
the best fit to the actual data, along with the la intrinsic width. The red lines show the equivalent from the shown simulated data set. 
Note how the CDJ model reproduces the sample properties best, followed by the CIH model. 



easily detectable. A systematic search for such shocks in 
clusters and galaxies is needed to assess whether all bub- 
bles are overpressured close to the center, but beyond the 
scope of this paper. However, for the overpressure model 
to work, we must also explain how all clusters overpres- 
sure their bubbles in a very similar fashion to stay in 
agreement with the correla tions discovered in this p aper. 
So far, onl y Centaurus A (iKraft et al.| [200l l2007f) and 
NGC 3801 (jCroston et al.ll2007D shOTTevidence of shocks 
surrounding their radio lobes. In addition, an initially 
overpressured bubble would have adjusted its bubble size 
on the order of a sound crossing time for the bubble, 
which is much shorter than the observed bubble ages. 

We also do emphasize that the current assessment of in- 
completeness effects may be oversimplified and influence 
these results. In order to properly address this important 
issue, an automated bubble detection tool is desperately 
needed. Incompleteness has a considerable impact on the 
energy budget associated with AGN feedback, and has to 
be taken into account in further studies of the subject. 



We also suggest that the current analysis be repeated 
with a larger and more homogenous sample, as the pre- 
dictions of the models will be much easier to distinguish 
from each other. 

We provide a framework to test not only currently 
available models, but that can easily be extended to fu- 
ture models of jet-lobe systems. We encourage modelers 
to find theoretical considerations predicting cavity sizes 
as a function of distance to cluster centers, or to ex- 
tract this information directly from their numerical sim- 
ulations. 



We thank Wise et al. for allowing us to reproduce their 
Hydra A figure, and Brian McNamara and Masanori 
Nakamura for stimulating discussion on the subject. This 
work was carried out in part under the auspices of the 
National Nuclear Security Administration of the Depart- 
ment of Energy (DOE) at Los Alamos National Labora- 
tory and supported by a LANL/LDRD program. 
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